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ABSTRACT 

We have studied the fluctuations of the soft (0.9-2 keV) X-ray background inten- 
sity for ^10 arcmin and ^2 arcmin beam sizes, using 80 high galactic latitude medium- 
deep images from the ROSAT position sensitive proportional counter (PSPC). These 
fluctuations are dominated (and well reproduced) by confusion noise produced by 
sources unresolved with the beam sizes we used. We find no evidence for any ex- 
cess fluctuations which could be attributed to source clustering. The 95 per cent 
confidence upper limits on excess fluctuations A/ C i us are: (A/ c i us //)io arcmin ^ 0.12, 
(A/ c i us //)2 arcmin < 0.07. We have checked the possibility that low surface brightness 
extended objects (like groups or clusters of galaxies) may have a significant contribu- 
tion to excess fluctuations, finding that they are not necessary to fit the distribution 
of fluctuations, and obtaining an upper limit on the surface density for this type of 
source. Standard Cold Dark Matter models would produce AI/I larger than the above 
limits for any value of the density of the Universe = 0.1 — 1 , unless the bias param- 
eter of the X-ray emitting matter is smaller than unity, or an important fraction of 
the sources of the soft X-ray background (~30 per cent) is at redshifts z > 1. Limits 
on the 2-10 keV excess fluctuations are also considered, showing that X-ray sources 
in that band have to be at redshifts z > 1 unless f2 > 0.4. Finally, if the spatial 
correlation function of the sources that produce these excess fluctuations is instead a 
power law, the density contrast Sp/p implied by the excess fluctuations reveals that 
the Universe is smooth and linear on scales of tens of Mpc, while it can be highly 
non-linear on scales ~ 1 Mpc. 

Key words: X-rays: general - X-rays: background - diffuse radiation - large-scale 
structure of Universe - methods: statistical 



1 INTRODUCTION 

Recent optical identification projects using ROSAT PSPC 
observations have resolved an important fraction of the ex- 
tragalactic X-ray Background (XRB) in the ~l-2 keV band 
into discrete sources, mostly Active Galactic Nuclei (AGN) 
and Narrow Emission Line Galaxies (NELGs, a mixed bag 
including Seyfert 2 galaxies, starburst galaxies and galaxies 
with HII regions) (Page et al. 1996a, Jones et al. 1995, Boyle 
et al. 1995, Boyle et al. 1994). 

Below 0.5 keV the fraction of the XRB that is extra- 
galactic is uncertain, with estimates ranging from about 10 
to 20 per cent (McCammon & Sanders 1990, Barber & War- 



wick 1994). The rest has a local origin, probably in a bubble 
of hot gas surrounding the Sun. Above 2 keV, only ~4 per 
cent of the XRB has been resolved. Ongoing identifications 
of serendipitous sources in ASCA images have increased this 
fraction to about 40 per cent (Inoue et al. 1996). 

Whatever the nature of the sources that produce the 
XRB, and independently of their identification, the inten- 
sity of the XRB received from different directions in the sky 
contains information on the angular distribution and clus- 
tering properties of such sources. The study of the distribu- 
tion of XRB intensities P(I) probes the source flux distri- 
bution (dN/dS or number of sources per sky area per unit 
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flux) down to fluxes S below the detection limit (Barcons 
et al. 1994, Hasinger et al. 1993). This technique is called 
P(D), (D — I — (Jxrb)) or fluctuation analysis and is most 
sensitive to fluxes in which there is about one source per 
'beam' (Scheuer 1974, Barcons 1992), the reason being that 
brighter sources contribute to the bright tail of the (skewed) 
distribution, while fainter and more numerous sources pro- 
duce negligibly small noise. However, if the counting noise 
is important, the technique is only sensitive to source fluxes 
equivalent to the photon counting noise level. 

The effect of source clustering is to decrease the effective 
number of sources per beam, hence broadening P(I) (Bar- 
cons 1992). This broadening can be related to the clustering 
properties of the sources that produce the XRB, which in 
turn are due to density fluctuations in the Universe Sp/p 
(Butcher et al. 1996, Barcons & Fabian 1988, Rees 1980). 

Instead of following the usual approach of using the 
deepest fields available to push our knowledge of dN/dS 
well below the present detection limits, in this work we have 
explored the clustering properties of X-ray sources by mea- 
suring or limiting the excess fluctuations they produce. Di- 
rect deep source counts have been performed over small sky 
areas, and they might be biased by large scale fluctuations 
in the source counts. The use of 80 widely scattered ROSAT 
fields allows a statistical study to be made (through P(I)), 
avoiding any such biases. 

The limits obtained on the excess fluctuations are then 
compared with the specific expectations from a Cold Dark 
Matter (CDM) model, to constrain the density of the Uni- 
verse (fi = 2qo) and/or the bias parameter of X-ray emitting 
matter with respect to the underlying matter distribution 
(6x). Assuming instead a power-law shape for the spatial 
correlation function of the source of the soft XRB, the up- 
per limits obtained on the excess fluctuations have been used 
to investigate 8p/p on different scales. 

In Section 2 we describe the data used in this work and 
the reduction process. A brief summary of P(I) analysis is 
given in Section 3, along with the dN/dS models used and 
the results of fitting the theoretical P(I) curves to the data. 

Section 4 is devoted to the development of the theoret- 
ical framework necessary to relate these excess fluctuations 
to CDM power spectra and 5p/p . The limits obtained on 
f2 and fex are also presented and discussed, as well as those 
obtained on 8p/p. In Section 5 we summarize our results. 

We have parametrized the Hubble constant as Ho = 
100/ikms -1 Mpc~\ with h = 0.5. The X-ray fluxes S will 
be given by default in the 0.5-2 keV range. 



2 THE DATA 

The data used in this work consist of 80 ROSAT PSPC 
pointings with exposure times longer than 8 ks at galactic 
latitudes higher than 20°. These same fields were used for 
the RIXOS survey (Mason et al., in preparation). In ad- 
dition, the RIXOS fields were chosen avoiding extended or 
very bright targets (e.g. clusters, nearby bright galaxies and 
bright stars). 

The Starlink software package ASTERIX was used for 
the data reduction. The data were screened for high particle 
background intervals (Plucinsky et al. 1993), bad aspect ra- 
tio solutions, and total accepted count rates deviating from 



the average of each observation. This procedure normally re- 
duced the nominal exposure time by 10 to 20 per cent. The 
remaining particle background was then calculated using the 
formulae in Plucinsky et al. (1993), and subtracted. 

The remaining counts in Pulse Height Analyzer (PHA) 
channels 92 to 201 (~ 0.7 — 2 keV) for each pointing were 
then binned to obtain images with a pixel size of 4.5 arc- 
sec. These images were then devignetted by dividing by the 
exposure maps provided by the standard EXSAS process- 
ing, after normalizing the maps to unity in the centre. We 
note however that the results given below are practically in- 
sensitive to whether the remaining particle background is 
subtracted or not, or on whether the vignetting has been 
corrected for or not. 

The range of channels used in this work was chosen 
to avoid local contributions to the XRB (such as the local 
bubble and Galactic diffuse emission, both thought to be im- 
portant only below ~1 keV), solar contamination (usually 
modelled as an oxygen line at about 0.5 keV, Snowden & 
Freyberg 1993) and absorption from neutral hydrogen (prac- 
tically absent above 1 keV) . An estimate of the possible so- 
lar contamination was obtained by extracting images just 
using night time observations (Snowden & Freyberg 1993). 
This reduced dramatically the total number of counts, hence 
worsening the statistics, without actually changing signifi- 
cantly the average count rate. We have, therefore, used both 
day and night time data. 

A circle of radius 5 arcmin around the target of each of 
the PSPC fields (generally at the centre) was excluded. This 
proved sufficient to exclude contributions from the targets 
down to the level in which their 'tails' would contribute less 
than 30 per cent of the local background per pixel in the 
two worst cases. In most of them this contribution was < 5 
per cent. 

Only one of the detected sources in the analyzed area in 
these fields is above the flux interval used in our calculations 
(see Section 3.1), excluding that field from our analysis does 
not affect any of our results, therefore we have used the 
80 fields including the detected sources within the regions 
explained below. 

Counts in each of the devignetted, particle-background 
subtracted, target-subtracted images were further grouped 
in two beam sizes: 

• An annulus of radii 5 and 10 arcmin centred on the 
pointing direction (the inner radius is due to the target sub- 
traction), giving a beam size of Q. c g — it x (10 2 — 5 2 )/3600 = 
0.06545 deg" 2 . The distribution, P(I), of the 80 XRB inten- 
sities obtained (in counts per second per 'beam') is shown 
in Fig. 1. 

• Eight circles of radius 2.5 arcmin with centres equally 
spaced in a circumference of radius 7.5 arcmin centred 
on the pointing direction, hence Qf ff = n x (2.5/60) 2 = 
0.00545 deg~ . We excluded two of these circular beams be- 
cause more than one third of their area was taken away by 
the target exclusion circle (that was slightly off centre). The 
remaining 638 values (80x8 —2), 7s, again in cts -1 beam -1 , 
give P(Is), as shown in Fig. 2. 

Both sets of intensities cover similar detector zones, but 
they sample different angular scales: 10 to 15 arcmin in the 
first case and <5 arcmin in the second. The maximum offaxis 
angle used (10 arcmin) ensures that the vignetting correction 
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Figure 1. Histogram of the distribution of the XRB intensities 
for the large beam (see text). Also shown as a solid continuous 
line is the best fit P(I) with K = 55 deg -2 , with no cluster 
contribution and the average A/ no ; se (see text). 
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Figure 2. Histogram of the distribution of XRB intensities for 
the small beam (whole dataset, see text). Also shown as a solid 
continuous line is the best fit P(Is) with K = 55 deg -2 , with no 
cluster contribution and the average A/snoisc (see text). 

is small (< 5 per cent) and that the effective area is also 
uniform over the detector region used. 

We found average values of the XRB intensity of 
0.49 ± 0.02 ct s _1 deg" 2 from the large beam sample and 
0.50 ± 0.03 cts -1 deg -2 from the small beam sample (both 
1 sigma confidence intervals). We adopt (7xrb) = 0.50 ± 
0.03 cts" 1 deg" 2 . 

The photon counting noise was estimated by the 
square root of the number of counts in each 'beam' (us- 
ing poisson statistics). We found AJ no i se = 0.0018 ± 
0.0005 cts" 1 beam" 1 and A7 S noi sc = 0.0005 ± 0.0002 
cts" 1 beam" 1 , in both cases we give 1 sigma confidence in- 
tervals. 

A conversion factor of lets" 1 (92-201) = 2.02 x 
10" 11 erg cm" 2 s" 1 (0.5-2 keV) was used throughout, ac- 
curate within ~5 per cent for power-law energy spectral 
indices a ~ 0.4 — 0.7, hydrogen column densities A<h ~ 



(0.5 — 20) x 10 20 cm" 2 and any combination of detector re- 
sponse matrix and effective area, thus covering the observed 
XRB spectrum (Gendreau et al. 1995, Branduardi-Raymont 
et al. 1994) and the galactic columns of the ROSAT obser- 
vations used (Mason et al. 1996). 

We therefore measure a total XRB intensity (including 
sources) of (7 X rb) = (3.3 ± 0.3) x 10" 8 erg cm" 2 s" 1 sr" 1 
(0.5-2 keV). This value is somewhat higher than previ- 
ous XRB intensity estimates, but still overlaps within ~2 
sigma with the value obtained by Branduardi-Raymont et 
al. (1994), for example. 



3 FLUCTUATION ANALYSIS 

3.1 Contribution from point sources 

In this work we have adopted the dN/ dS shape and param- 
eters from Barcons et al. (1994): 

dN k ( sy- 

= £ (£) s>Sb 

with S B = 2.2 x 10" 14 erg cm" 2 s"\ 7 d = 1.8, 7 U = 2.5 and 
K =55 deg" 2 . 

The results given below do not change if we use the 
slightly different parameters from Branduardi-Raymont et 
al. (1994) or Hasinger et al. (1993), which is hardly surpris- 
ing considering that they are all mutually consistent, have 
been obtained with ROSAT data and sample similar or over- 
lapping flux ranges. This also means that no biases have 
been introduced in the determination of the source counts 
in those surveys by large scale source number fluctuations. 

The dN/dS parameters given above are appropriate for 
S between 0.07 and 50x 10" 14 erg cm" 2 s" 1 . At the level of 
one source per beam, the P(I) curve is going to be sensitive 
down to fluxes S~4x 10" 14 erg cm" 2 s" 1 and the P(h) 
down to S ~ 0.5 x 10" 14 erg cm" 2 s" 1 . The sensitivity limit 
of our analysis is thus S ~ 10" 14 erg cm" 2 s" 1 . The width 
of the P(I) is mainly due to this 'confusion noise' rather 
than to photon counting noise. Although we are integrating 
the dN/dS between zero and infinity in our calculations, the 
practicalities of using a Fast Fourier Transform algorithm to 
calculate the P(I) effectively reduced this interval to (0.02 — 
40) x 10" 14 ergcm" 2 s" 1 . 

Given a dN/dS and a beam profile, the shape of P(I) 
can be predicted (see Barcons 1992 and references therein). 
The counting noise is generally taken into account by con- 
volving P(I) with a gaussian of width A7 no i se . We have also 
followed this approach, taking as A7 no i se the average values 
given above, and checking the influence of the dispersion 
around those values by using the 1 sigma upper and lower 
limits as well (see below). 

The beam functions are taken as two circular step func- 
tions: one with an outer radius of 10 arcmin and an inner 
radius of 5 arcmin (for 7), and another with just an outer 
radius of 2.5 arcmin (for 7s). The sizes are much larger than 
the Point Spread Function (PSF) of the XRT/PSPC com- 
bination (Hasinger et al. 1992), making the convolution of 
the PSF and the step functions indistinguishable from the 
simple step functions in practice. 
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Any width in excess of that expected from the source 
flux distribution and the poisson counting noise is called ex- 
cess variance, and is usually modelled by convolving P(I) 
with a gaussian of width AJ c i u8 . We assume that the ex- 
cess fluctuations arise from clustering of sources, and per- 
haps some contribution from extended sources like clusters 
of galaxies (see below). If any other unknown systematic ef- 
fect contributes to the excess fluctuations, the results given 
below would just be upper limits to A7 c i us really due to clus- 
tering, and any consequences of the results given here would 
be strengthened. 

The model P(I) is then (see also Eq. 22 in Barcons 
1992): 



P(I) 



duj e 



01 exp {- 



-uj 2 AI 2 



c /2-w 2 A/ 2 lus /2} 



"eff- 



x exp jfi cff / dSdN/dS[exp(2TTiioS/n ctt ) - l] j (1) 
The same expression is valid for 7s replacing Q B g with 
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Figure 3. Contours of AL values (1, 2 and 3 sigma) in (ii",AJ G i us ) 
space, for the large beam, with no cluster contribution and the 
mean A7" no ; se (see text) 



3.2 Contribution from extended sources 

We have considered the X-ray emitting clusters reported 
by Rosati et al. (1995). With the above parametrization, a 
set of parameters that follow their dN/dS in the flux range 
(1 - 40) x 10- 14 ergcm -2 s- 1 is: 7d = 7 „ = 1.962, K c \ = 
9.784 deg~ 2 and SW = 1CT 14 erg cm" 2 s" 1 . 

The properties of these clusters have been taken from 
the study of poor groups of galaxies by Mulchaey et al. 
(1996). We have assumed the temperature of the hot gas (re- 
sponsible for the detected X-ray emission) to be kT ~ 1 keV 
and a King emission profile with a cluster core size of 
7?corc = 15 arcmin (changing the size to 7 arcmin does not 
affect the results given below) . For a nearby group (like those 
in Mulchaey et al. 1996) with z ~ 0.02 this corresponds to 
a core size of ~ 0.4 Mpc (or ~ 0.2 for 7 arcmin). 

The conversion factor for clusters, assuming a thermal 
bremstrahlung spectrum with the above temperature and 
absorption by neutral hydrogen with A^h = 10 20 cm -2 , is 
lets -1 (92-201) = 1.64 x 10" 11 erg cm" 2 s" 1 (0.5-2 keV). 

The cluster contribution to P(I) is modelled by con- 
volving it with the P(I) due to the clusters only, i.e. by 
adding another term to the exponent in braces in Eq. |l| 



P(I) 



duj e 



x exp < 



x exp 



nwI exp {-o; 2 AI 2 oisc /2 - u 2 AJ 2 lus /2} 
dSdN/dS[exp(2iriu)S/Q. eS ) - l] j 



2nJdrrJdS(dN/dS) c i [exp(2mu}SG cl (r)) - l] 1(2) 

where G c i(r) is the convolution of a King profile with the 
step functions described above. 

Only clusters with fluxes S > 10~ 14 erg cm" 2 s" 1 (the 
sensitivity limit of the Rosati et al. sample) have been used 
to calculate the P(T). Again, the results given below do not 
change if we decrease this limit by a decade, because of the 
flatness of the clusters source counts. 

The angular size of the clusters of galaxies considered 
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Figure 4. Contours of AL values (1, 2 and 3 sigma) in 
(if, ATsdus) space, for the small beam, with no cluster contri- 
bution, the mean Aig no j se and the whole dataset (see text) 



here implies that, if one of them is present in a given ROSAT 
pointing, the eight small beams will be affected. This intro- 
duces a correlation between them and complicates the error 
estimates on A/sdus- A way around this problem is to select 
one of the eight beams for each ROSAT pointing at random, 
and just use those 80 values of Is- This allows us to estimate 
the significance of the cluster contribution (at the price of 
sacrificing sensitivity). Should this contribution prove to be 
negligible, the whole dataset can be used, applying Eq. [j] 
instead of Eq. ^| 

3.3 Fitting process and results 

A Maximum Likelihood fitting method was adopted, x 2 was 
not adequate because the number of fitting points for P(I) 
was too small to make a significant number of bins with a 
reasonable number of points in each one of them (enough 
for gaussian statistics to be valid). 

A further ingredient (apart from dN/dS, (dN/dS) c i, 
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A/noise and A7 c i ua ) is necessary to fit P(I): the total in- 
tensity of the sources in the dN/dS used to calculate P(I) 
{{P/dN/ds) is always smaller than the mean observed inten- 
sity. The missing sources are not important for the shape 
of P{I), because there are so many of them and they are so 
faint that they contribute a negligibly small gaussian noise to 
it (already taken into account with A7 no i se ). Their absence 
makes the 'peak' of the model P(I) to be at an intensity 
smaller than that of the peak of the observed distribution, 
so an overall shift of the distribution is necessary to compare 
the observed and modelled P(I). 

An additional intensity A7 is added to each 7 to shift 
them to higher values and is allowed to vary until a best 
fit is obtained (keeping the rest of the parameters fixed). 
It is then discarded as a non-interesting parameter and the 
fitting proceeds with a different set of parameters. The best 
fit A7 is in fact only very weakly dependent on the rest of 
the parameters. The value of A7 can however be predicted 
from the dN/dS and (7xrb), and its final best fit value is 
not expected to be very different from this predicted value. 

For each set of fitting parameters, we have defined the 
likelihood function as 




^clus ( cts S 1 beam ) 

Figure 5. One dimensional AL profiles extracted from the con- 
tours in Fig. 3 (see text) as a function of AJ c i us . The solid line 
corresponds to the mean A/ no ; se , and the dashed and dotted lines 
correspond to adding and subtracting la from it, respectively. 



L(A7 clus ,70 = -2^mP(7i) 



+ 



(7xrb) — A7 — (I)dN/ds\ 



A(7 X rb) 



(3) 



where the first term is the usual definition (and P(I) is as 
defined in Eq. [l] or []), and the second term makes added 
intensities far from their expected values less likely, weighted 
for each beam size with the error in the estimate of the XRB 
intensity, A(7xrb), given above. With this definition AL is 
distributed as Ax 2 - 

The first fit is performed fixing all the dN/dS parame- 
ters to the values given above and leaving A7 c i us as the only 
free parameter. The best fit values are shown in Table 1 (for 
the large beam) and Table 2 (for the small beam) . The effect 
of the uncertainty on A7 no ise has been assessed by fixing it 
to its mean value and the 1 sigma upper and lower limits, 
and performing the fit for each of these three values. The 
results are indicated in Tables I and 2 (rows with both the 
K and K c \ columns labelled 'fixed'), with the first row of 
each group of three corresponding to the mean, and the sec- 
ond and the third line to the I sigma upper limit and lower 
limit, respectively. At the 2 sigma confidence level, only up- 
per limits are obtained: A7 c i us < 0.004 ct s -1 beam -1 and 
A7 S cius < 0.0005 cts -1 beam" 1 (or A7 c i us /(7 X rb) <12 per 
cent and A7s c 1u S /(7xrb) <19 per cent). 

The dN/dS normalization, K, and A7 c i us are coupled 
to some extent: large normalizations increase the 'intrinsic' 
P(7) width, thus reducing the amount of excess variance 
needed. We have done a second set of fits in two dimensions, 
with both K and A7 c i us as free parameters. The results are 
shown in Tables I and 2. AL contours are plotted in Fig. 
3 (large beam) and 4 (small beam), for the case of no clus- 
ter contribution (and the whole dataset, see below) and the 
average values of A7 no lse and A7s n oisc, respectively. 

It is possible to obtain confidence intervals on A7 c i us 
taking into account its coupling with K by finding the 
minimum AL value as a function of K for every A7 c i us , 
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Figure 6. One dimensional AL profiles extracted from the con- 
tours in Fig. 4 (see text) as a function of A7g c i us . The solid line 
corresponds to the mean A/g no ; sc , and the dashed and dotted 
lines correspond to adding and subtracting la from it, respec- 
tively. 



and then considering them as a one dimensional AL pro- 
file for A7 c i us (Lampton, Margon and Bowyer 1976). This 
has been done for the results plotted in Fig. 3 and 4, 
and it is shown in Fig. 5 and 6, respectively, as well as 
in Tables 1 and 2. As in the one dimensional 
the 2 sigma confidence level, only upper limits are ob- 
tained: A7ci us < 0.004 cts -1 beam" 1 and A7 Sc ius < 0.0006- 
0.0008 cts - beam - (or A7 c i us /(7xrb) <12 per cent and 
A7 S c1u S /(7 X rb) <22-30 per cent). 

All the above fits have been repeated without any clus- 
ter contribution, and the results also included in Tables 1 
and 2. It is clear that adding the clusters does not sig- 
nificantly reduce the L values, nor does it reduce the ex- 
cess variance. We obtained a quantitative assessment of the 
significance of this contribution using the standard F-test 
(Bevington 1969). This assesses the relative improvement in 
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Table 1. Results of the fit to P(I). 



A T 

J noise 


A T 

^clus 

I K^li D UCdilll J 


2cr 

(ct s 


upper limit 


T 

L 


A 

/j„„-2\ 
K ae & ) 


K cl 
(Ape,- 2 ) 


2cr upper limit 
v ae S ) 


Mean 


0.0017 




0.0037 


-556.3 


55.0 fixed 


9.8 fixed 


- 


Mean+lcr 


0.0008 




0.0034 


-556.3 






- 


Mean— lcr 


0.0021 




0.0039 


-556.3 


" 


" 


- 


Mean 


0.0006 




0.0038 


-556.8 


55.0 fixed 


25.8 


62.2 


Mean+lcr 


0.0000 




0.0035 


-556.6 




19.3 


55.8 


Mean— lcr 


0.0013 




0.0039 


-556.8 


" 


26.0 


67.2 


Mean 


0.0001 




0.0038 


-557.0 


67.3 


9.8 fixed 


- 


Mean+lcr 


0.0000 




0.0037 


-556.7 


61.9 






Mean— la 


0.0003 




0.0040 


-557.1 


71.1 






Mean 


0.0022 




0.0041 


-555.6 


55.0 fixed 


0.0 fixed 




Mean+lcr 


0.0016 




0.0038 


-555.6 








Mean— lcr 


0.0025 




0.0043 


-555.6 








Mean 


0.0001 




0.0038 


-557.1 


74.1 


0.0 fixed 




Mean+lcr 


0.0003 




0.0037 


-556.5 


62.9 






Mean— lcr 


0.0001 




0.0039 


-557.2 


78.3 







Table 2. Results of the fit to P(I S ). 
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0.0007 


-862.0 


49.5 


9.8 fixed 






Mean+lcr 


0.0000 




0.0006 


-861.8 


31.0 








Mean-lcr 


0.0003 




0.0008 


-861.9 


55.0 






80 


Mean 


0.0000 




0.0004 


-862.1 


55.0 fixed 


0.0 fixed 






Mean+lcr 


0.0000 




0.0004 


-858.6 










Mean-lcr 


0.0003 




0.0006 


-862.2 








80 


Mean 


0.0000 




0.0007 


-862.2 


52.5 


0.0 fixed 






Mean+lcr 


0.0000 




0.0006 


-862.0 


34.1 








Mean-lcr 


0.0003 




0.0008 


-862.2 


55.0 






638 


Mean 


0.0000 




0.0002 


-6820 


55.0 fixed 


0.0 fixed 






Mean+lcr 


0.0000 




0.0002 


-6797 










Mean-lcr 


0.0003 




0.0004 


-6820 








638 


Mean 


0.0000 




0.0002 


-6823 


62.6 


0.0 fixed 






Mean+lcr 


0.0000 




0.0002 


-6800 


47.3 








Mean-lcr 


0.0000 




0.0002 


-6836 


76.5 







X (or L) on the addition of a new free fitting parameter 
(.Kci); in our case that means comparing the values of L in 
the second and fourth groups of rows in Tables ^ and ^| 
The F-test reveals that the addition of does not im- 
prove significantly the fits, to a confidence of 96 per cent 
for the large beam, and the best fit for the small beam ac- 
tually corresponds to K c \ = 0. We can then conclude that 
their contribution to the P(I) width is negligible and ignore 
it. This allows us to use the full 638 values of Is, reducing 
considerably the 2 sigma upper limit in the excess variance: 
A/sdus < 0.0002 cts -1 beam -1 (or AJ S c1 US /(/xrb) <7 per 
cent). 

Confidence regions on K c \ can be obtained from the 



AL contours in the (AJ c i us , K c {) space with the method de- 
scribed above. Only upper limits are obtained at 2 sigma 
level, and are given in Tables ^ and ^. Rosati et al. warn that 
their value is only a lower limit to the real surface density of 
clusters (or extended X-ray sources). Our results show that 
down to 10 -15 - 10 -14 erg cm s , the surface density of 
clusters is not larger than 3 to 6 times the value obtained 
by Rosati et al. 

Soltan et al. (1996) found an important contribution 
(~30 per cent) to the angular correlation function of the 
soft XRB from extended haloes around Abell clusters of 
galaxies on scales > 1 degree. Since we are exploring much 
smaller angular scales and the opposite (low flux) end of the 
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dN/dS distribution of the X-ray emitting clusters, there is 
no contradiction between our finding that extended sources 
(clusters) do not contribute significantly to the excess fluc- 
tuations and the results of Soltan et al. (1996). 

The upper limits on the excess fluctuations obtained in 
this section (namely, A7 c i us <12 per cent and A/sdus <7 
per cent, with a 2 sigma confidence level), will be used in 
Section 5 to constrain the values of the density parameter 
of Q and 6x using the expressions derived in Section 4. 



4 INHOMOGENEITIES IN THE MASS 
DISTRIBUTION OF THE UNIVERSE 

4.1 Relation of excess fluctuations to the power 
spectrum 

It is easy to realize that (A7 c i U s/(7xrb)) 2 is the value of 
the autocorrelation function of the XRB at zero-lag. We 
can then use the expressions in Appendix A of Barcons & 
Fabian (1988) and Eqs. 2 and 4 of Carrera et al. (1991) 
to relate the limits found on the excess fluctuations to the 
clustering properties of the sources of the XRB. 

In our case, the beam shape is a two sided step function, 
with a value of 1 between ri and r 2 and outside, where 
ri = 5 arcmin and r 2 = 10 arcmin for the large beam, and 
n = and r 2 = 2.5 arcmin for the small beam. Its two 
dimensional Fourier transform is 



4.2 Luminosity functions and modelling 



G(q) = (r 2 Ji(r 2 q) - r 1 J 1 (r 1 q))/q 



(4) 



where q is the magnitude of the two dimensional Fourier 
space vector, and Ji(x) is the Bessel function of order 1. 

Solving Eq. 2 and 4 in Carrera et al. (1991) for (7xrb), 
we arrive at 



1 

P 



A/ C l us 

(7xrb) 



4V2tt 



x w 1 dz ( l + z )~ a ( l + 2 ^r 1/2 j 2 ^)/d 2 A (z) 



d 2 qG 2 {q)i(q/d A {z)) 



dz{l + z)~ t '{l + 2qoz)- 1/2 j(z) 



(5) 



dA(z) being the angular distance and j(z) the K-corrected 
volume emissivity (emitted power per unit volume) of the 
sources that produce the excess variance and contribute a 
fraction / to the XRB. £(fc) is the three dimensional Fourier 
transform of the spatial correlation function £(r), and k is 
the magnitude of the three dimensional Fourier space vector. 
Following Peebles (1980), £(fc) is also the power spectrum, 
multiplied by (2/tt) 3 ^ 2 , due to the different definitions of the 
Fourier transform used here and in Peebles (1980). 

Eq. ^ allows the calculation of A7 c i us /(7xrb) for a par- 
ticular 6x and a power spectrum model, which in turn would 
depend on fi (see below). By comparing these predictions 
with the upper limits obtained above, constraints can be 
placed on those cosmological parameters. In the next sec- 
tion we present the luminosity function we have used to 
calculate j(z). 



At the flux levels at which our P(7) anaysis is sensitive 
(S ~ 10~~ 14 erg cm -2 s _1 ) the dominant type of X-ray 
sources found in ROSAT surveys are AGN although with 
an increasingly important fraction of NELGs (McHardy et 
al., in preparation, Mason et al., in preparation, Boyle et al. 
1995, Carballo et al. 1995, Boyle et al. 1994). 

The X-ray Luminosity Function (XLF, number of 
sources per unit volume and unit luminosity) of AGN has 
been very well studied recently with ROSAT at those fluxes 
(Page et al. 1996a, Boyle et al. 1994). It has been found to 
be well represented by a broken power law. Within a pure lu- 
minosity evolution model, the AGN luminosities have a fast 
positive evolution up to z ~ 1.5 — 2. At that redshift the 
evolution slows down, or even stops and becomes negative. 

We have obtained the emissivity j(z) in Eq. |E] by in- 
tegrating the best fit XLF models of Page et al. (1996a), 
since AGN are the main contributors to the XRB over the 
flux range studied. Indeed, AGN are about ~50 per cent of 
the sources at the fluxes we are dealing with, and we have 
to consider the redshift evolution of the volume emissivity 
from other sources. The evolution of NELGs, the other type 
of source with a sizeable contribution and likely to be clus- 
tered, is somewhat different (Page et al. 1996b, Boyle et al. 
1995). Their rate of evolution is lower than that of AGN, 
and they are concentrated at low z (< 0.6). 

We have adopted the best power-law model with cut 
off evolution and go = 0.5 with a conversion factor of 1.8 
(between ROSAT and Einstein fluxes) from Page et al. 
(1996), but making qo half the value of f2 investigated in 
each case. The other best fit models produce very similar 
A7 c i us /(7xrb) values. The K correction has been calculated 
using a ~ 1, as observed for AGN, the dominant type of 
sources in our flux range (Mittaz et al., in preparation, Al- 
maini et al. 1996, Ciliegi et al. 1996, Romero-Colmenero et 
al. 1996, Vikhlinin et al. 1995). 

We have also considered the results on 2-10 keV 
Ginga excess fluctuations from Butcher et al. (1996): 
(A7 c iu S /(7xrb)) < 0.038 (2 sigma). In this case the red- 
shift dependence of the emissivity j(z) of the sources is not 
known and we have adopted a very simple model for their 
redshift distribution: j(z) oc (l + z) i+p . p — corresponds to 
no evolution of the emissivity in comoving coordinates. For a 
simple power law luminosity function, p ~ 3 implies a lumi- 
nosity evolution similar to that found in the soft band. We 
have approximated the Ginga collimator shape by a gaus- 
sian of dispersion 7s ~ 0.8 deg, and used an energy index of 
a = 0.7 as observed for AGN in that band. 

Making / = 1 in Eq. |B| is equivalent to assuming that the 
sources whose clustering produces the excess fluctuations we 
are studying produce all the XRB. We know that only 50-60 
per cent of the XRB is produced by sources with fluxes larger 
than ~ 10 -14 erg cm -2 s _1 (our sensitivity limit). However, 
since in Eq. ^ the absolute normalization of the XLF cancels 
out, just having more sources with the same evolution would 
not affect our theoretical A7/7. We have also checked that 
extending the integrals in redshift in Eq. ^| to z = 5 instead 
of 2 = 3 (our default value) does not affect our results. If, as 
discussed above, the NELGs are proved to make an impor- 
tant contribution to the XRB, but with a different evolution 
(more concentrated at lower z), the resulting density fluc- 
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Table 3. A7" cIus /(/ X RB> from Cold Dark Matter for ROSAT. 
Linear power spectrum from Peacock & Dodds (1994) 



n 


Uxrb) 


A/ Clus 
Uxrb) 


0.1 


0.107 


0.098 


0.2 


0.098 


0.086 


0.4 


0.097 


0.081 


0.(5 


0.101 


0.080 


0.8 


0.106 


0.081 


1.0 


0.111 


0.083 


Upper limits 


0.072 


0.119 



tuations produced by these sources would be larger, hence 
strengthening our results. 

A similar argument can be used for the Ginga upper 
limits. 

4.3 CDM Models 

Cold Dark Matter models present a picture of the Universe 
in which the smallest structures (galaxies) form first and, by 
merging, form larger structures (Peacock & Dodds 1994). 
Even if the basic assumptions have not been thoroughly 
tested, the CDM scenario provides useful calculation tools 
and expressions to analyze the evolution of the Universe. 

This is the case for the power spectrum of density fluc- 
tuations P(k). A number of useful parametrizations that fit 
some of the available angular and spatial clustering data are 
found in the literature (Peacock & Dodds 1994 -hereafter 
PD-, Efstathiou, Bond and White 1992 -hereafter EBW-, 
Bardeen et al. 1986). 

We have used the shape of the linear power spectrum 
of PD: 

^1^1+^) I 2 
v ' 47rfc 3 \ gk J 

x | [l + ak + (bk) 2 + (ckf + {dkf] ~ 1/4 j 2 (6) 

where a = (3.89/P) h' 1 Mpc, b = (14.1/r) hT 1 Mpc, 
c = (5.46/r) h' 1 Mpc, d = (6.71/r) ft" 1 Mpc, g = 
(2.34/P) ft -1 Mpc, and T is a shape parameter that can be 
changed, both to make Eq. ^ fit several different observa- 
tions, and to reflect the behaviour of different CDM and 
Mixed Dark Matter models. Following PD, we have chosen 

r = fi/iexp(™2fi B ) (7) 

which is equivalent to that presented by EBW for zero 
baryonic density Qb = 0, but also includes an empirical 
dependence in Qb, making high baryonic content models 
mimic low CDM density. The power spectrum parametriza- 
tion with the shapes and parameters from EBW is similar 
(for a power spectrum index n = 1). 

PD also give a dependency of the normalization of P(k) 
with O: as = 0.75f2 -0 ' 15 , erg being the rms density con- 
trast when averaged over spheres of radius 8 ft -1 Mpc. We 
have adopted this normalization dependence on Q. For each 
value of we have calculated as from Eq. ^, rescaling its 
normalization to give the value of <xg(f2) given above. This 
normalized P(k) is then used to calculate AI /I. 



Table 4. AI c \ us /{Ixrb) from Cold Dark Matter for Ginga. Lin- 
ear power spectrum from Peacock & Dodds (1994) 



Zmax 


P 


= 


V 


= 3 


1 


3 


1 


3 


Q 




A-f c l us 


/Oxrb) 




0.1 


0.219 


0.222 


0.066 


0.113 


0.2 


0.174 


0.176 


0.051 


0.081 


0.4 


0.142 


0.143 


0.040 


0.059 


0.6 


0.129 


0.129 


0.036 


0.050 


0.8 


0.121 


0.122 


0.033 


0.045 


1.0 


0.116 


0.117 


0.032 


0.043 


Upper limit 


0.038 


0.038 


0.038 


0.038 



Standard primordial nucleosynthesis and abundances 
observations constrain Qb ~ 0.05 (Olive & Steigman, 1995), 
and we have assumed this value. Since Ob only appears in 
an exponent and is small in any case, changing it by ±0.01 
(its observational confidence interval) does not change the 
results given below. 

X-ray sources are possibly more clustered than the un- 
derlying matter, and therefore the A7 c i us /(/xrb) obtained 
from CDM has to be multiplied by the bias parameter 6x- 
A value 6 X ~ (3.4 ± 0.8)fi ' 6 //' has been found for nearby 
bright X-ray sources, where /' is the fraction of the gravi- 
tational acceleration on the Local Group contributed by the 
z < 0.015 region (/' ~ 0.5, Miyaji 1994). 

With all the above assumptions, our CDM 
A7 c iua/(/xR B ) only has two free parameters: Q and 6x- We 
have sampled f2 between 0.1 and 1, and assumed 6x = 1. 
Different values of £1 change the shape of the CDM power 
spectrum (through the shape parameter P), while the effect 
of 6x is just multiplicative. 

The above CDM power spectrum shape is constant in 
comoving coordinates. Its evolution with redshift is obtained 
by multiplying its normalization by a factor D 2 (z) that is 
proportional to (1 + z)~ 2 for Q, = 1, and has a more com- 
plicated dependence with redshift for smaller values of Q. 
(Peebles 1980). 

We present in Table ^ A/ c i us /(/xrb) obtained for the 
beam sizes and shapes used here (small and large beam), for 
several different values of il in the above range, bx = 1 and 
the PD power spectrum given in Eq. ^| (the power spectrum 
of EBW produces similar results). AI /I produced by CDM 
exceeds our small beam upper limits (A7s c 1us/(/xrb) < 
0.07,) for any value of fi. The power spectrum of the spatial 
distribution of the X-ray emitting matter is not compatible 
with CDM. 

We have also used the non-linear scaling of the power 
spectrum proposed by PD. This only increases the excess 
fluctuations from CDM (by about 50 per cent for the small 
beam, the more stringent limit), hence worsening the mis- 
match. A faster clustering evolution does not therefore help 
reconcile CDM with the excess fluctuations upper limits. 

If either &x < 1 (i.e., the X-ray sources are less clustered 
than the underlying mass distribution) or / < 1 (i.e., the 
sources considered in our XLF do not produce the whole of 
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the XRB), CDM models would be consistent with our excess 
fluctuations upper limits, provided that / x b x ^ 0.7. 

We have already shown that the X-ray sources more 
clearly associated with peaks on the matter distribution 
(clusters) are not relevant for the excess fluctuations. How- 
ever, AGN and NELGs have been shown to be important 
contributors to the soft X-ray background (50 to 60 per cent 
of it has been resolved into these types of sources), and both 
populations seem to cluster in the same comoving scales as 
'normal' galaxies do (see Boyle & Mo 1993 for a study of the 
clustering of X-ray AGN, Shanks & Boyle 1994). Values of 
6x between 1 and 8 would be obtained from the results of 
Miyaji (1994), with Q varying in the above range. A value 
of the bias parameter bx < 1 is therefore very unlikely. 

As discussed at the end of Section 4.2, the absolute 
normalization of the emissivity of the sources that produce 
the excess fluctuations j(z) cancels out. We also commented 
that, if an important fraction of those sources were dis- 
tributed at smaller redshifts than the population considered 
in the XLF used here, the calculated excess fluctuations pro- 
duced would increase. About 95 per cent of the excess fluctu- 
ations from CDM are produced at z < 1; sources at higher 
redshift do not contribute significantly to the excess fluc- 
tuations. From this it follows that a possibility of getting 
/ ~ 0.7 to reconcile CDM and XRB fluctuations would be 
to place the unresolved part of the sources of the XRB at 
z > 1. 

AIcIus/{Ixrb) calculated for a Ginga beam size and 
a power law emissivity evolution are given in Table ^ for 
p — 0,3 and two different values of the maximum redshift 
of integration z m&K = 1,3. The minimum redshift was set at 
0.05, changing it to 0.1 did not change the results signifi- 
cantly. For a comoving evolution p — 0, the upper limits are 
exceeded at all values of Q, and the maximum fraction con- 
tributed by 2 < 1 sources is / < 0.3. A positive evolution is 
in principle more plausible, in line with the soft XLF results 
quoted above. For p = 3, about 50 per cent of the XRB in- 
tensity has to come from z > 1 to reconcile the upper limits 
with the CDM excess fluctuations. Alternatively, most (70- 
90 per cent) of the XRB sources could be nearby, but then 
the density of the Universe cannot be low (fi > 0.2 — 0.4). 

Similar upper limits on / were obtained from studies 
of the angular correlation function of the XRB both above 
and below 2 keV and in angular scales between 1 arcmin 
and several degrees (see e.g. Carrera et al. 1991, Soltan & 
Hasinger 1994). However, the alternative possibility in those 
studies of a rapid evolution of the source clustering would 
not be consistent with our data, as discussed above. 



4.4 Limits on the density contrast from excess 
fluctuations 

In this section, we will investigate the density contrast of 
matter in the Universe (Sp/p) implied by the upper limits 
obtained on the excess fluctuations. Instead of using a CDM 
power spectrum, we assume that the sources of the XRB 
have a spatial correlation function £(r) = (r/ro)~ 1,s with 
a comoving evolution. By performing its Fourier transform 
and substituting in Eq. [| we can translate the limits on 
A/ c i us /(7xrb) to limits on the spatial correlation length ro. 
The density contrast Sp/p in a window W(r) is given 

by 



o 
o 



Q. 



°0.1 



Ginga . 

_d i i 

1 10 100 

Comoving scale (Mpc) 



1000 



Figure 7. Sp/p limits at different scales (see text): the solid line 
is the ROSAT upper limit from the large beam, the dashed line 
is the small beam upper limit, the dotted line is the Ginga upper 
limit. 



(Sp/p) 2 



d 3 k£(k)W 2 (k) 



(8) 



where W'(k) is the Fourier transform of the window function, 
that we have taken here to be a sphere of radius R. For 
this window function and a power law correlation function, 
Sp/p is also a power law on R: (Sp/p) 2 oc (ro/R) 1 ' 8 . We 
can therefore use the limits on ro obtained from Eq. |B| (with 
the emissivities discussed in Section 4.2) to constrain Sp/p. 
The resulting upper limits on Sp/p versus R are plotted 
in Fig. 7, using both our limits on the excess fluctuations 
from ROSAT, and Butcher et al. (1996) results from Ginga 
(assuming p — 3 that gives the more conservative upper 
limits). We have used Q — 0.1 in Fig. 7. If instead we use 
= 1, the limits are 10-20 per cent smaller. 

Given the size of the different beams used, this analysis 
is going to be sensitive to different sampling radii. We have 
estimated the relevant ranges by using a typical angular dis- 
tance for each beam size involved (~ 3' for our ROSAT small 
beam, ~ 12' for our large beam, and ~ 1° for Ginga) and 
calculating the maximum and minimum separations it cor- 
responds for the redshift range considered (z ~ 0.05 — 3 for 
ROSAT, and z ~ 0.1 — 3 for Ginga). As we can see in Fig. 
7, the smaller angular scale results are sensitive to spatial 
distances of the order of 1 Mpc, while the larger ones are 
sensitive to a few tens of Mpc. 

At the larger scales sampled here the Universe is quite 
homogeneous (Sp/p < 1), while below a few Mpc there is 
space for strong density fluctuations (Sp/p > 1), that would 
reveal a highly non-linear growth of structure. 



5 SUMMARY 

Our fluctuation analysis of 80 ROSAT fields has allowed us 
to constrain the excess fluctuations on ~10 arcmin angular- 
scales to be A/cius < 0.004 cts -1 beam -1 and on ~2 arcmin 
A7s c iu S < 0.0002 cts" 1 beam" 1 (or AJ c1us /(J X r B } <17 per 
cent and A7sc1us/(^xrb} <7 per cent), both with 2 sigma 
confidence levels. 
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The source counts found in medium and deep surveys 
in empty fields reproduce well the fluctuations of the XRB 
around bright targets (most of which are nearby galaxies 
of different types). Since there is no need for any excess 
fluctuations, we conclude that faint X-ray sources are not 
associated to local astronomical objects. 

A contribution from extended objects with low surface 
brightness (like groups or clusters of galaxies) is not required 
to fit the observed distribution of intensities. The surface 
density of these objects is shown to be <3 to 6 times the ob- 
served value, limiting the fraction of low surface brightness 
sources missed by present surveys. 

The upper limits on A7 c i us /(7xrb) obtained here (and 
others from the literature) have been compared with CDM 
theoretical models to extract constraints on the density pa- 
rameter of the Universe f2 and the bias parameter of X- 
ray emitting sources with respect to the underlying mat- 
ter distribution &x- Unless bx ~ 0.7 (which is unlikely), 
the only possibility for reconciling our results with CDM 
would be that the remaining unresolved sources of the soft 
XRB (contributing 30 per cent of it) are at z > 1, and have 
suffered a cosmological evolution different from the other 
known sources of the soft XRB (AGN and NELGs). Simi- 
larly, sources that produce about 50 per cent of the 2-1 keV 
XRB have to be at z > 1; this fraction could be larger if 

n > 0.4. 

In a different approach, a power-law shape is assumed 
instead for the spatial correlation function of the XRB 
sources, constraining the density contrast to be < 1 on 
scales of tens of Mpc and < 10 — 100 around one Mpc. This 
indicates that the Universe is very homogeneous at larger 
scales, but inhomogeneities might be present and common 
at smaller scales, as observed in surveys of the nearby Uni- 
verse. 
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